16.451 Lecture 23: A potential model for nuclei Nov 28, 2005 



From last class: 

a simple square well potential V(r) 

• We have a complicated system of A 
nucleons. 

• About half of them are protons, so a 
repulsive (+ve energy) term has to 
be added to the square well to 
account for this (~ few MeV) 

How to connect this model to something 
observable? 


Independent particle model: 

• Assume independent particle motion in some average nuclear potential V(r) as shown. 

• Then we can fill the eigenstates of the potential to maximum occupancy to form a 
nucleus, as is done with electrons in atoms (to 1 st order...) 


Connection to average nuclear properties: 2 

• The binding energy of each nucleon, in our model, is a few MeV. 

• The potential energy of a bound nucleon is negative, by ~ 0.3% of its rest mass energy, 

which therefore has to show up as a decrease in its mass. 

• For A nucleons, the total binding energy is: 



B= Z B ' = M 


mass of nucleus, M 


The average binding energy per nucleon, B/A, can be determined from mass data 
and used to refine a model for V(r); it ranges systematically from about 1-9 MeV 
as a function of mass number for the stable isotopes. 


Reference: F&H ch. 16 
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Atomic Mass Units: 
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• By convention, we set the mass of the carbon-12 atom as a standard. 

• Denote atomic masses with a "script" M, measured in atomic mass units, U 

M( 12 C) = 12.0000000000 U (exact!) 1 U = 931.494 MeV (expt.) 


Calculation for carbon-12: 

m p = 938.2 MeV ' 
m n = 939.6 MeV y 
m e = 0.511 MeV , 


6xV m, = 11,269.8 MeV 


12 U = 11,178.0 MeV 

I 


B( 1 6 2 C) = ^ nij-M = 91.8 MeV 


Binding energy per nucleon in 12 C: B/A = 7.8 MeV; 
Contrast to the deuteron 2 H: B/A = 1.1 MeV 


The famous Binding Energy per Nucleon curve: 4 
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Implications of the B/A curve: 


Greater binding energy implies lower mass, greater stability. 

Energy is released when configurations of nucleons change to populate the 

larger B/A region -> nuclear energy generation, e.g. 


Fusion reactions at small A release substantial energy because the B/A curve rises 
steeply at small A: 


2 H+ 2 H => 3 He + n + 3.27 MeV V 
=> 3 H + p + 4.03 MeV > 
2 H + 3 H => 4 He + n + 17.6 MeV ; „ 


Binding energy 
of products is greater 
than the sum of the 
binding energies 
of the initial species. 


Fission reactions at large A release energy because the products have greater 
binding energy per nucleon than the initial species: 

235 U + n => 100 X + 133 Y + 3n + 200 MeV I 

v V J 

distribution of final products 


A semi-empirical model for nuclear binding energies: 


1. Volume and Surface terms: 

First consider a 1-dimensional row of nucleons with interaction energy per pair = s 
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total: A 

\ 1 


each has 2 neighbors correction 

for the ends 


B r, 

— = 2s - 

A 

A 
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Approximately constant, with end 
effects relatively smaller at large A. 


By analogy, for a 3-d nucleus, there should be both volume and surface terms 
with the opposite sign, the surface nucleons having less binding energy: 


73 /i A 2/3 . B /i — 1/3 

-D — CZy^iT. Qg/\. — s ~r — Q-y CZ § 
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2. Coulomb term: 
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for a uniform sphere, 


-'Coul 


r q(r) dq 
J 4 ns 0 r 


total charge + Ze 


3 (Zef 
5 4ns 0 R 


This effect increases the total energy and so decreases the binding energy. 
Simple model: AB = - a c Z 2 A _1/3 


But this is not quite right, because in a sense it includes the Coulomb self energy 
of a single proton by accounting for the integral from 0 to r p ~ 0.8 fm. The nucleus 
has fuzzy edges anyway, so we will have to fit the coeff icient a c to mass data. 


Solution: let AB scale as the number of proton pairs and include a term: 

AB = -a c Z(Z -1) A _1/3 => A5- = - a c Z (Z -1) A“ 4/3 


3. Symmetry Term: 


So far, our formula doesn't account for the tendency for light nuclei to have Z = N. 
The nuclear binding energy ultimately results from filling allowed energy levels in a 
potential well V(r). The most efficient way to fill these levels is with Z = N: 


Simplest model: identical nucleons as a Fermi gas, i.e. noninteracting spin- i 
particles in a box. Two can occupy each energy level. The level spacing ~ 1/A. 

A mismatch between Z and N costs an energy price of AE at fixed A as shown. 


energy 
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r J here.. 





AB = — a A {Z — N) 2 A -1 = -a A (A- 2Z) 2 A -1 
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4. Pairing Term: 


Finally, recall from slide 1 that for the case of even A, there are 177 stable 
nuclei with Z and N both even, and only 6 with Z and N both odd. Why? 

-> Conf igurations for which protons and neutrons separately can form pairs 
must be much more stable. All the even-even cases have J n = 0 + , implying 
that neutrons and protons have lower energy when paired to total angular 
momentum zero. 


Solution: add an empirical pairing term to the binding energy formula: 


II 
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with +1 for even-even, 0 for even-odd, and -1 for odd-odd 

Full expression: 

B(Z,A ) = a v A-a s A 2/3 ~ a c Z (Z -1) A“ 1/3 - a A ( A-2Z ) 2 A -1 + 5 


Fitting of coeff icients to data: 1° 


B(Z, A) = a v A-a s A 2l3 -a c Z (Z-l) A _1/3 — a A (A-2Z) 2 A -1 + 8 
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Semi-Empirical Mass (binding energy) Formula (SEMF) implications: 12 


B(Z,A) = a v A-a s A 2/3 -a c Z(Z-l)A~ m - a A (A-2Z) 2 A -1 + 8 


Stable nuclei have the maximum B for a given A; for constant mass number, 
B is quadratic in Z -> “mass parabolas", e.g.: 
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Beyond the SEMF: "Magic Numbers" and the Shell Model 


We already noted that there 
were some marked deviations 
from the SEMF curve at small 
mass number, e.g. A = 4. 

On an enlarged scale, a systematic 
pattern of deviations occurs, with 
maxima in B occurring for certain 
“magic" values of N and Z, given by: 

N/Z = 2, 8, 20, 28, 50, 82, 126 

These values of neutron and proton 
number are anomalously stable 
with respect to the average - the 
pattern must therefore ref lect 
something important about the 
average nuclear potential V(r) that 
the neutrons and protons are 
bound in.... 

(NB, the most stable nucleus of all is 
56 Fe, which has Z = 28, N = 28, 

"double magic" ...) 



Other evidence for "magic numbers" 2-n and 2-p separation energies: 14 



• Energy required to remove a pair of 
neutrons or protons from a given nucleus is 
referred to as S 2n or S 2p 

• like the ionization energy for atoms, but 
the pairing force is so strong in nuclei that 
systematics are more easily seen comparing 
nuclei that differ by 2 nucleons 

• the same pattern of “magic numbers" 
appears - large separation energies 
correspond to particularly stable nuclei: 

N/Z = 2, 8, 20, 28, 50, 82, 126 .... 
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A periodic table of nuclei? 
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Systematics are reminiscent of the periodic structure of atoms, which results from 
filling independent single-particle electron states with electrons in the most eff icient 
way consistent with the Pauli principle, but the magic numbers are different: 


“Magic numbers" for atoms: 

Z = 2, 10, 18, 30, 36, 48, 54, 70, 80, 86 .... 



z 


Figure 5,1 Atomic radius (top) and ionization energy (bottom) of the elements, 
The smooth variations in these properties correspond to the gradua! filling of an 
atomic shed, and the sudden jumps show transitions to the next shell. 


The Nuclear Shell Model 
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Self-consistent approximation: assume the quantum state of the i th nucleon can be found 
by solving a Schrodinger equation for its interaction via an average nuclear potential V N (r) 
due to the other (A-l) nucleons: 


-ft n 

— v 2 + V N (r) 

2ju 


Wntmifi) =E nl y nlm {r i ) 


Assume a spherically symmetric potential V N( r); 
then the eigenstates have definite orbital angular 
momentum, and the standard radial and angular 
momentum quantum numbers (n,/,m)as indicated. 

(Justification: measured quadrupole moments of 
nuclei are relatively small, at least near the 
“magic numbers" that we are interested in explaining; 
midway between the last two magic numbers, ie 
around Z or N = 70, 100, the picture changes, and 
we will have to use a different approach, but at 
least for the lighter nuclei this assumption should 
be reasonable.) 
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Shell Model continued... 
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If we choose the right potential f unction V N (r), then the wave f unction for the whole 
nucleus can be written as a product of the single particle wave functions for all A 
nucleons, or at least schematically: 

^ Nucleus ( F ) = n=l Vnlmin) 


oversimplification here... actually, it has to be written as an 
antisymmetrized product wavefunction since the nucleons are 
identical Fermions - the procedure is well-documented in 
advanced textbooks in any case! 

With total angular momentum given by: 

A 

J = X ^ ^ = ^ + ’ ( s = 7 ) 

/=! 


And parity: 



Always + for an even 
number of nucleons... 


What to use for V N (r)? - three candidate potential functions: 
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Square we) I | 

r< ft: V{r)~ - V g 
: l/fri=0 


Harmonic oscillator 

m = 


Saxon-Wood s 


V{r)= 


1+BX p 


Advantage: easy to write down 
Disadvantages: 

numerical solutions only 
edges unrealistically sharp 


Advantage: easy to write down and 
can be solved analytically 
Disadvantage: potential should 
not go to infinity, have to cut off 
the function at some finite r and 
adjust parameters to fit data. 


Advantage: same shape as measured 
charge density distributions of 
nuclei, smooth edge makes sense 

Disadvantage: 

numerical solution needed 
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Comparison: Harmonic Oscillator versus Woods-Saxon solutions: 19 


• since both potentials are spherically 
symmetric, the only difference is in 
the radial dependence of the wave 
functions 

• amazingly, when parameters are 
adjusted to make the average 
potential the same, as shown in the 
top panel, there is remarkably 
little difference in the radial 
probability densities for these 
two potential energy f unctions! 

• this being the case, the simplicity 
of the harmonic oscillator potential 
means that it is strongly pref erred 
as a model for nuclei 



Evidence that this works: 
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z 

electric charge density, measured via electron scattering: P$) = e Vi (0 1 2 

i = 1 


charge density difference 
between 205 TI and 206 Pb is 
proportional to the square 
of the wave f unction for 
the extra proton in 206 Pb, 
i.e. we can actually measure 
the square of the wave 
function for a single proton 
in a complex nucleus this 
way! 


But, we still have a problem 
explaining the magic numbers 
- next class! 
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